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Features of the nuclear isothermal incompressibility n and adiabatic incompressibility kq are 
investigated. The calculations are done at zero and finite temperatures and non zero entropy and for 
several equations of state. It is shown that kq decreases with increasing entropy while the isothermal 
K increases with increasing T. A duality is found between the adiabatic kq and the T = isothermal 
K. Our isothermal results are compared with a recent lattice Monte Carlo calculation done at finite 
T. The necessity of including correlations is shown if k is to have a peak with increasing T as seen 
in the Monte Carlo calculations. A peak in k is linked to attractive scattering correlations in two 
- ; - nucleons channel in the virial expansion in our approach which are Pauli blocked at low T. 
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, The behavior of nuclear systems at moderately high temperature and density is of current interest for several 
reasons. Such studies are important for understanding features of current medium energy collisions 1], for future 
bJO RIA experiments, and for nuclear astrophysics as in supernovae explosions. The incompressibility and its associated 
^ ' equation of state of pressure versus density and temperature is important in understanding flow produced in nuclear 
I collisions as reviewed in Ref.0 and, in general, the thermodynamic properties of fermionic system under these 
conditions. This paper focuses on an important quantity for understanding properties of these systems which is the 
nuclear incompressibility. While the nuclear incompressibility at zero temperature has been studied for an extended 
. period Jjt is only relatively recently that its temperature dependence has been of concern. A recent investigation 

^ ' of Miiller etal [3| using quantum Monte Carlo techniques with nuclei on a lattice studied this temperature dependence 
, and showed some interesting features about it. Namely, they have shown that the incompressibility coefflcient can 
become very large at temperatures around 14 MeV. Large increases in the specific heat at constant volume was also 
seen in their quantum Monte Carlo approach. This increase in Cy is also seen in Ref. 7] using a totally different 
approach based on recursive methods to obtain the finite temperature partition function of hadronic matter. In this 
later approach inhomogeneities in the nuclear density are incorporated in the model in the form of clusters. When a 
homogeneous system breaks into clusters at relatively low density, cost functions associated with large surface energies 
^ arc present and are responsible for the dramatic increase in Cy with T. 

, In this paper we will study the behavior of the isothermal incompressibility with T and properties of the equation 
' I ■ of state (EOS) at higher than normal density. As a baseline, we will begin with a mean field discussion of its behavior 
O [ with T to see how large the incompressibility can become and if reaches the large values and associated peak found 
2 in the Monte-Carlo result of Ref.Q. 
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First, we define a quantity k, the incompressibility coefficient, as 

■ , d\E/A) ^ d\E/A) 

This quantity is evaluated at the saturation density po where E / A has a minimum. The giant monopole resonance 

energy is then i?o = \J'^^^- If the temperature is kept constant in the above derivatives we have the isothermal 

incompressibility k, and if the entropy is held fixed, the result is the adiabatic incompressibility k,q. The k and kq are 
equal at T = only. The quantity k defined above is not the isothermal compressibility defined in thermal physics as 

K = --{^){dV/dP)T with T held fixed and here P is the pressure. Since P = -dF/dV, we have K = (i) (l/^) ■ 

At T = 0, F = — TS = E, and thus k = 9/ (pqK). This reciprocal connection between K and k is no longer true at 
finite T. Besides the isothermal compressibility, an adiabatic compressibility Kq = —{y){dV/dP)s can be obtained 
by keeping the entropy constant. The reciprocal is related to kq, the adiabatic incompressibility as kq = 9/{pqKq). 
Since the natural variables for energy are entropy and volume from dE — TdS — PdV variations of the energy with 
V at constant S bear a similar relation to variations of the Helmholtz free energy with V at constant T where 
dF = —SdT — PdV . Thus the adiabatic incompressibility of Eq.Q can go to zero for a Skyrme interaction as we shall 
see. The minimum point (also maximum point) in the energy occurs at zero pressure since dE/dV at constant entropy 
is —P. Therefore E at constant S has the same maximum and minimum points with variations in V or R or density 
as F at constant T since both derivatives are —P which is set to 0. The behavior of the adiabatic incompressibility 
is linked to a phase change. As we shall see, the behavior of the isothermal incompressibility may be associated with 
the appearance of a strongly corelated system at high density and temperature. 

Our mean field discussion is based on a Skyrme interaction. To keep the discussion simple, we consider uncharged 
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symmetric nuclear matter with no surface energy terms. The Skyrme interaction energy then is 

[//A = -aop + a„p^+". (2) 

The ao term gives a medium range attraction while the Oq, term is a short range repulsion. 

At T = 0, the kinetic energy Ek/A = {1)Ef{p) with the Fermi energy Ep = ^(^f-p)^^''- The coefficients oq 
and Oct are fixed to give a binding energy per particle Eb/A = 16MeV at density p = po = 0.15/fm"^, which gives 
^0 = o-oPo = 37 + 23/a and ba = AqPo^" = 23/a in MeV. The incompressibility coefficient k at T = is then 

K = -2Ek/A + 9(1 + a)aaapl^°' = 165 + 207a. (3) 

For a ~ 1/3, k — 234MeV. Smaller values of a lead to softer equations of state and lower k. In the limit a — > 0, 
logarithmic terms appear in Eq.© coming from the presence of a factor (x/a)(l — a;") — a;log(a;). The x = 
p/Pq = (Rq/R)^ with Rl = A/{po4:Tt/3). The a ^ limit is the softest EOS allowed by Eq.Q, and this limit gives 
from Eq.Q a value of k = 165 MeV. A stiff EOS has a = 1 and k = 372 MeV. Recent calculations done at T = 
0,0, [l3 have a value of k = 210 ~ 270 MeV and suggest a value of a = 1/3 in a Skyrme type approach. A larger 
range of values of k, from 211 ~ 350 MeV, were reported in Rcf. llj. Because of these uncertainties in k, from more 
realistic forces, we will present results for various values of a, from a = to a = 1. 

At non-zero T <C Ep, the kinetic energy is Ek/A — i^)Ep + {^)T^/Ep 12]. The energy per particle is 

E/A = 21xV3 + (^) 1^ - box + (4) 

The value of x that minimizes E/A is Xm and satisfies the equation: 14(a;m'^ — x^^") — &o(a^m — ^In") — i^)'^'^ / ■ 
Then the k — k{T) is given by 

k{T) = -A2xl(^ + 0.705TVa;2/3 + 9a{l + a)bc,xl;t°' (5) 

At T = 2.5, 5, and 7.5MeV, and for a = 1/3, the values of k are 242, 265 and 302 MeV, respectively. The corresponding 
values of Xm are: 1.011, 1.043, 1.091. When T is replaced with entropy per particle S/A then this T dependent term 
becomes {Ep/7r^){S/Ay = {35/n^){S/A)^x'^/^ since S = {tt^/2)TA/Ep at low T. This S/A term can simply be 
added to the first term on the right side of Eq. ^ since both have the same x^^"^ dependence. 

If the corrections to the nuclear matter incompressibility at T = from finite temperature terms are small, then 
these corrections can be obtained by using the following method. Let Eo{R) be the nuclear matter energy per 
particle EOS and which has a minimum at Rq and an incompressibility kq. If we add to this a term E^{R), so that 
E{R) = Eo{R) + Ex{R) then the minimum shifts to a new point Rm = i?o + ^Rx- The new minimum and k can be 
found by making a Taylor expansions around of Rq. The new n is 
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The various quantities are evaluated at i?o and each ' represents one derivative wrt R. Corrections to k involving the 
skewness Sk — R^{Eq{R) / A)'" or third derivative of the energy were pointed out in Ref.^jQ. Ellis et al used the 
correlation between compression modulus and skewness coefficient to examin the implications in a relativistic Hartree 
Fock approximation where the E^ is the Coulomb interaction. The above expression is a modified version of their 
result. Eq.® gives an expression for the incompressibility at T = 0. This will be kq = ^0(0;) in Eq.©. The skewness 
is Sk{a) = -3(509 + 828a + 207a^) in MeV. Comparing Sk with k of Eq.Q we see that the ratio of Sk/n^ is of the 
order of 10 and somewhat insensitive to a. 

At low T, tdking E:,{R,T) = OmUT'^R^ /A^/^ , we obtain the K{a ^ 0) = 165 + 1.16T2, ^(a = 1/3) = 234+1.32r2 
and K{a = 1) = 373+1. 61T^. Thus we see that the first term is very sensitive to a but the finite temperature correction 
is somewhat insensitive to a. At fixed entropy, the second derivative of E{R)/A has a very different behavior than 
at fixed T. Namely, it decreases with S/A. This can easily be seen by noting that Ex{R, S) = A.%i&{S/ A^ A^/^ / R^ 
compared to E^{R, T) = Q.^bllT'^R^ / A^/^ . We have the following final results (in MeV): K(a -> 0) = 165-30(5/71)2, 
K(a = 1/3) = 234 - 38(S'/A)2 and K{a = 1) = 373 - b2,{S/Af. 

At higher T, the nearly degenerate Fermi gas kinetic energy term is replaced by a virial expansion in pX^ , where 

A = is the quantum wavelength. Namely, Ek/A = (3/2)r(l + Y.n Cn(pA3/4)") with coefficients arising from 

antisymmetrization that are: ci = 1/2^/^ = 0.177, C2 = (1/8 - 2/3^/^) = -3.3 x 10"^ C3 = 1.11 x 10"", • • • [Tlll^. 
Since the c„'s become small rapidly, we will keep terms up to C2. Then k is given by: 

k{T) = -27r(AVo/4)'c2x2„ + 9a(l + a)5„xi,+" (7) 
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The Xm is again the minimum of E/A, but now evaluated with the new kinetic energy. The Xm is affected by both 
the Ci and C2 terms at temperatures where Ci dominants. A limiting value of k can be obtained by taking T very 

large where Ci term leads a minimum Xm given by = ^ij^"^^, (^1 — i^f^^'^) • this high T limit k is given by 

the second term on the right side of Eq. and is 

niT ^ lOMeV) = K^at (l - ^T^^^) (8) 

and goes to its saturation value Ksat = 9q;6o(6o/(1 + a)&Q)^^" with a T dependence of l/VT. We note that the sign 
of ci determines whether it approaches from above or below. For purely antisymmetric correlations ci is positive 
because of the statistical repulsion of fermions. If Ci becomes negative as will be discussed below it would approach 
from above. At infinite T for a — 1, k, = 704MeV with Xm = 1.304 or a minimum density pm ~ 1.304po and for 
a = 1/3 K = 468MeV with a;,„ = 1.53 or = X.hZp^. In the limit a 0, k = 380 MeV with x,„ = e"/23 ^ ^.84 
or pra — 1.84/90. Thcse are the limiting values for k and Xm- The mean field results with antisymmetrization effects 
only can be compared with the Monte Carlo results of Ref.[a| which show a very sharp increase in the k with T until 
a temperature of ^ 14 MeV is reached. Then the incompressibility sharply decreases. A peak in k, of about 1500 
MeV is present at a saturation density of 2po = 0.3/fm'^. From T = 20 MeV to higher T's the incompressibility is flat 
or saturates at a value of about 250 MeV. It should also be noted that the saturating value of k obtained by Monte 
Carlo are lower than ours even for the softest EOS with a limit. The difference is due to finite size surface effect. 
The lattice calculations are done for small systems. An explanation of the maximum is qualitatively given in terms 
of clusters that form which then repel each other through next nearest neighbor interactions which is repulsive. In 
our mean field model the incompressibility increases with T and saturates at a value which is about twice its value 
at T = 0. The saturation density in our model also increases with T and is somewhat smaller than the Monte Carlo 
results of 2po- 

Before discussing how a peak in k may arise in our approach, we briefly investigate the case of constant entropy 
in the ideal gas limit using the Sackur- Tetrode law 0: SjA = 5/2 — ln(A'^p/4). This law connects T to p or F as 
T = Gsp^l'^- Here Cs = I (mc^)\ exp[(2/3)(S'/A) - 5/3]. The resulting E(R)IA has a structure similar to the 

result for a degenerate Fermi gas since both have a p^l^ dependence for the kinetic energy term but with different 
coeflficients. This feature and a similar result at lower T suggests a duality in the energy per particle EOS at constant 
entropy and its associated kq and the T = EOS and its associated constant T k. We also note a parallel between 
_F as a function of T and V and i? as a function of S and V. 

We now turn to the issue of clusters or more precisely correlations at moderately high T and high p > pQ. We 
study the corrections to the ideal gas law using the virial expansion P — pT (l + ci(pA'^/4) + C2(pA'^/4)^ + ■-). If 
antisymmetry effects are the only corrections, the coefficients can be calculated by following a procedure in Eef . fwi Il7l | 
and are the coefficients already given before Eq.Q. This procedure is an extension of our fragmentation model by 
simply noting that a cycle of length k is analogous to a cluster of size k with a weight function Xk- The grand 
canonical partition function is \ogZ = Y^Xke^^^ and the mean number of cycles is (n^) = xue^^^ The pressure 
is PV/T = \ogZ = 'Y^XkC^^^ . A constraint exists, {A) = 'Y^k{nk) = A which determines the fugacity z = e^'^ in a 
power series in A by inverting the series. Then we arrive at 

PV , —X2 ,9 4x9 — 2x1X3 ,0 
- A + — -A^ ^ — ' A^ 



T xi 
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-20x^ + 18x1X2X3 - 3x^X4 ^ 



Substituting Xk = (— l)*^^^(y/A"^)/fc^/^ gives the desired power series in {A/V)X^/4: for fermions [Ulil- The factor of 
4 is spin and isospin degeneracy. The same procedure applies for bosons with Xk — {V/ X^)/k^^^ . For fragmentation in 
the Boltzmann limit the Xk = xH^k) ^int{k) = xi[k) ^^^^^ where Fk is the internal free energy of a cluster of size k and 
A(fc) = X/k^^^ . The effect of antisymmetry for odd k clusters and symmetry for even k clusters can be included. The 
grand canonical ensemble represents a system of fermions (odd cluster sizes) obeying FD statistics and bosons (even 
cluster sizes) obeying BE statistics. The constraint of chemical equilibrium Mi = ku^OY Oh = zpLp + n/i„ is imposed 
which determines the fugacity from the constraint. In the x^ model of Ref.jl^, ll^j H3 this amounts to having 
various terms in Xk that represent both cycles and clusters. For example xi = p-, X2 = —-^j^-^ + 2^/^p-Zi„t(2)/4 
and X3 = '^nj^ + 3'^/^-p-Zi„t(3)/4. X4 will have terms from the antisymmetry of monomers from cycles of length 4, 
from symmetrization of dimers and from clusters of size 4. Once the x^'s are given the canonical partition function 
can be generated by a recurrence relation |l8l l22j | . A factor 1 /4 appears from spin- isospin degeneracy which has been 
included. The internal partiton Zint{k) = 9{E])e^^'^''^ + i Ej,t / ^^e-f^^dE. The sum is over 
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bound state Ej which have degeneracy g{Ej) and Sj^t is the phase shift in channel of spin J and iospin T. These 
phase shifts include effects from both attractive and repulsive interactions. Using nucleon-nucleon phase shifts the 
continuum contributions reduces the bound state contribution by about 50% for moderate temperature {T ^ 20 
MeV) and less for low temperatures because of the Boltzmann weight factor in the integral. At infinite T Zint 
since the continuum exactly cancels that bound states by Levinson's theorem 23] . As an initial example for Zi„t(2)2'^/^ 
we will consider (l/2)(3/4)2'^/^el'^^l/"^ to see how it compare with 1/2^/^; 1/2 is for the continuum reduction, 3 is 
for the spin degeneracy of the ground state of the deutron. This choice has a value 1.06 neglecting Eb- To reduce 
1.06 to 1/2^/^ we would need a reduction factor or 1/6. Also other spin isospin channels increase Zint- Thus ci can 
easily become minus. For a negative ci, the n is above its saturating value and approaches it from above as T~^/^. 
At low T, K is below its saturating value and initially increases as because the Fermi sea blocks excitations. This 
behavior automatically implies a peak in k. 

Since we are at higher than normal density, we do not think of real physical clusters, but now view the presence of 
k — 2 objects as correlated pairs of fermions in various isospin and spin scattering channel which are allowed since 
the Pauli blocking has been removed. A deutron like structure may appear as a resonance and the e'^^'/^ — > e^^"/^ 
with Er the resonance energy. Higher order fc = 3, 4, 5, ... terms represent higher order correlations of fermions in 
various J, T channels. For example k = Z can represent J = 1/2, T = 1/2 correlations. 

In this paper we investigated the behavior of the infinite nuclear matter incompressibility at finite temperature 
and entropy uing a mean field theory and also considering the role of correlations. Various forms of the EOS are 
studied using a Skyrme parametrization. Both the isothermal (constant temperature) and adiabatic (constant entropy) 
incompressibilities are found to be sensitive to the choice of the Skyrme repulsive parameter a which gives the power 
of the density involved in the repulsive term. These two incompressibilities have very different behaviors. The 
isothermal incompressibility increases with T initially as until a saturation value is reached while the adiabatic 
incompressibility decreases with increasing entropy and eventually goes to zero. The isothermal incompressibility 
approaches its saturation value as — with the minus sign reflects an approach from below the saturating value 
with increasing T. This behavior arises from the statistical asymmetric repulsive correlations that represent the Pauli 
exclusion principle. The adiabatic incompressibility is shown to arise from an equation of state (EOS) or energy per 
particle that has a structure that is similar to a T = Fermi gas. 

Our isothermal incompressibility results are compared with a recent Monte Carlo calculation which gave a much 
stronger temperature dependence and much larger value of the isothermal incompressibility. The Monte Carlo result 
has a peak behavior in k with a peak of k = ISOOMeV at T '--^ 14MeV. The saturation point is at twice normal density 
in the Monte Carlo result which is somewhat larger than the mean field saturation density. 

We then discussed how a peak can appear in the isothermal incompressibility by looking at coefficents in the virial 
expansion, and in particular, we investigated the first coefficient called ci here. The approach of k to its saturation 
value was shown to be related to the sign of ci , with ci positive having an approach from below and ci negative having 
an approach to the saturating value from above it. The role of attractive correlation between nucleons was studied and 
shown to most likely be strong enough in various spin isospin channels to overcome the statistical repulsion term at 
high T . The Pauli blocking of scattering is reduced at high T and such correlations are present even a density p > pQ. 
A two nucleon correlation of paired fermions in a high density, but also high temperature, medium can account for 
the existence of a peak. The structure of the peak is related to the presence of strongly correlated fermions in triplets, 
and higher order correlations. 
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